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More often than not, models of flavor symmetry rely on the use of nonrenormalizable operators 
(in the guise of flavons) to accomplish the phenomenologically successful tribimaximal mixing of 
neutrinos. We show instead how a simple renormalizable two-parameter neutrino mass model of 
tribimaximal mixing can be constructed with the non-Abelian discrete symmetry T-j and the gauging 
of B — L. This is also achieved without the addition of auxiliary symmetries and particles present 
in almost all other proposals. Most importantly, it is verifiable at the Large Hadron Collider. 
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In 2001, the non-Abelian discrete symmetry A4 was 
shown for the first time |l[ to allow for the seemingly 
incompatible pattern that charged-lepton masses are all 
very different and yet a symmetry exists to predict the 
neutrino mixing matrix without knowing the individual 
neutrino masses. In 2004, it was shown for the first 
time that A4 could also predict neutrino tribimaximal 
mixing with sin 2 29 atm — 1 and tan 2 8 so i = 1/2. Since 
early 2005, when the solar angle in neutrino oscillations 
was revised by SNO Q to tan 2 8 sol = 0.45 ± 0.05, this 
idea became widely accepted and the use of non-Abelian 
discrete symmetries [J for understanding flavor has ap- 
peared in very many publications [f|. The two earliest 
papers @, @] after the SNO revision in 2005 both used 
A4 and suggested two different two-parameter neutrino 
mass matrices, whereas the original proposal [2] of 2004 
had three parameters. 

If the 3x3 Majorana neutrino mass matrix M. v is 
rotated by the Cabibbo-Wolfenstein unitary matrix [H, [s| 



(1) 



where u> = exp(27ri/3) = —1/2 + i\/3/2, then the tribi- 
maximal form was shown Q to be 
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The two examples mentioned above are then [6j b = 
and 1 71 



UMvU 1 = 



-d 2 /a s 
ad 
da, 
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However, these forms are only obtained at the expense of 
additional auxiliary symmetries and particles, and with 
the use of nonrenormalizable operators @ . 

On the other hand, it has been shown recently [l(| that 
A4 alone is sufficient to obtain b = 0, if the alternative 



A4 lepton assignments of Ref. [ll| are used instead of the 
original proposal of Ref. [l[ and that neutrinos become 
massive through Higgs triplets [l2[ in a renormalizable 
model. Here we show how Eq. (3) may be obtained by 
the canonical seesaw mechanism for neutrino mass, using 
the non-Abelian discrete symmetry T7 [TJ| and gauging 
B — L [14], without the addition of auxiliary symmetries 
and particles or the use of nonrenormalizable operators. 

Most importantly, our proposal is verifiable at the 
Large Hadron Collider (LHC). The predicted Z' gauge 
boson will decay into scalars which support the T 7 sym- 
metry. Their subsequent decays into charged leptons will 
then reveal the predicted T7 flavor structure used in ob- 
taining neutrino tribimaximal mixing. 

Since there are three families, non-Abelian discrete 
symmetries with irreducible three-dimensional represen- 
tations are of special interest. The smallest group with 
a real 3 representation is A4 which has 12 elements. The 
smallest group with a complex 3 representation is T7 
which has 21 elements. The group A(27) [l5| is slightly 
bigger (27 elements) and also has a complex 3 represen- 
tation. They are all subgroups of SU(3). The various 
irreducible representations of the three groups are 



A(27) 



It (i = 1,2,3), 3; (4) 
U (i = 1,2,3), 3, 3; (5) 
1, (» = 1,2, 3, 4, 5, 6, 7, 8, 9), 3, 3. (6) 



Their crucial differences are in the following group mul- 
tiplications 



Ai 
T 7 



3x3 = ^l t 



(7) 



3x3 = 3 + 3 + 3, 3x3 = 3 + 3 + 3, 
3x3 = ^1^ + 3 + 3; (8) 

i 

A(27) : 3x3 = 3 + 3 + 3, 3x3 = 3 + 3 + 3, 

3x3 = ^1,. (9) 

i 

We will show that our T7 model assignments cannot be 
replaced by either those of A4 or A (27). 
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The finite group Tj is generated by two noncommuting 
3x3 matrices: 



(10) 



where p = exp(2iir/7), so that a 7 = 1, b 3 — 1, and 
ab = ba 4 . Let 3 = (x±, x 2 , x 3 ), and 3 = [x\, x 2 , x 3 ), 
then their possible multiplications form the following 
3 representations: {x 3 y 3l xiyi, x 2 y 2 ), (x 2 yi, x 3 y 2 , Xty 3 ), 
(x 2 ?/3±S3272,X3?7i±xiy3,xiy2±i22/i), and the following 
3 representations: (x 3 y 3 , Siyi, S2&2), (#12/2, x 2 y 3 , x 3 yi), 
(x2y 3 ±x 3 2/2,x 3 yi±xi2/ 3 ,xi?/2±a;22/i). The combinations 
xiyi +aj fe_1 X2y 2 + ui 2k ~ 2 x 3 y 3 form the representations 1^. 
for k = 1, 2, 3 respectively. 

Under T 7 , let L 8 = (i/,Z)j - 3, Z£ ~ L., i = 1,2,3, 
$ 4 = (0+, 0°); — 3, which means that $; = (0°, ~ 
3. The Yukawa couplings Lil^k generate the charged- 



lepton mass matrix 




u> 2 f 2 v 2 

U>f 2 V 3 
1 




(11) 



if V\ = v 2 — v 3 — v/\/3, as in the original A4 proposal [![. 

Let ~ 3, then the Yukawa couplings LiV^k are 
allowed, with 




(12) 



for v\ = v 2 = v 3 — t)/V3 which is already assumed for 
to/. Note that $ and l> have B - L = 0. 

Now add the neutral Higgs singlets Xi ~ 3 and 77^ ~ 3, 
both with B — L = —2. Then there are two Yukawa in- 
variants: vfVjXk and vfVjTjk (which has to be symmetric 
in Note that x* ~ 3 is not the same as to ~ 3 be- 
cause they have different B — L. This means that both 
B — L and the complexity of the 3 and 3 representations 
in Tj are required for this scenario. The heavy Majorana 
mass matrix for v c is then 



M 




(13) 







where A = hu\ — hu 2 = hu 3 and B = h'u' 2l u[ = ti 
have been assumed, i.e. Xi breaks in the (1,1,1) direction, 
whereas to breaks in the (0,1,0) direction. This is the 
Z 3 — Z 2 misalignment also used in A4 models. 



The seesaw neutrino mass matrix is now 
-moM' 1 - 



in 




,(14) 



which has only two parameters and is identical to Eq. (3). 
Detailed numerical analysis of this form was already done 
in Ref. Q. Here we achieve the same result without the 
auxiliary Z4 x Z 3 symmetry and extra particles assumed 
there. The key is that 3 x 3 x 3 is an invariant in TV, but 
not in A(27), whereas A4 cannot distinguish this from 
3x3x3 which yield two other invariants in T7. 

To realize the misalignment of (\) ~ (1, 1, 1) and (77) ~ 
(0, 1, 0), we need to choose the soft breaking terms in the 
Higgs potential consistent with these different residual 
symmetries [l6j]. However, the quartic terms XiXjVkVi 
have several T7 invariants, and most of them will destroy 
this pattern. To maintain the desired misalignment, this 
model has to be supersymmetrized. 

Consider x ~ 3 and 77 ~ 3 as supcrfields with B — L = 
-2. Add x' ~ 3 and rj' ~ 3 with B - L = 2. Then the 
superpotential contains the terms 

w = f^i^Xk + fijk^jVk + m xX ix'i + mnViV'i, (15) 
from which the F terms of the Higgs potential are 



Vf = \m xXt \ 2 + \m nVl \ 2 + \f^ kV f^ + m xX ' k \ 2 



(16) 



whereas the D terms from U(1)b-l are 

Vd = 2g 2 B _ L \ X *Xi + vtVi ~ xTx'i ~ €vT- (17) 



With the addition of bilinear soft terms x*Xii Xi x\i 
XiX'i + H.c., 77*773, 772*77^, 773772 + H.c, 77*771 + 773*773, 
Vi*Vi+V 3 *V 3 , mv'i+mv' 3 + H.c., and (xi + X2 + X3W2 + 
(xi + X 2 + X 3 ) r l2 + H.c, which preserve U{1)b-l as they 
must, T7 is broken with the desired pattern. 

Flavor-changing leptonic interactions through Higgs 
exchange are present in this model, but they are sup- 
pressed by lepton masses, as in the original A4 pro- 
posal [l|. The set of three Higgs doublets $i transform- 
ing as 3 under T7 is rotated by U of Eq. (1) to form 
mass eigenstates 4>o,i,2 ~ l 7 w,w 2 under the residual Z 3 , 
where 4>q is identified as the one Higgs doublet (with 
((/>q) = v) of the Standard Model, with Yukawa couplings 
v~ 1 [m e eLen + m^JiLPR + to t tltr], which is of course 
flavor-conserving. The "flavor-changing" interactions of 
0i,2 are then given by 

Ant = 7j 1 \m T L^ j^Tpi + m^LeLfiR + m e L TL eR]4>i + 
v~ x [m T L eL T R + m^LrLHH + m^L^^R^ + ff.c.(18) 

However, if the neutrino sector is ignored, a lepton 
flavor triality (Z 3 symmetry) (l7j exists here, under 
which e,/_t, r ~ l,w 2 ,w, implying thus the decays t + — > 



3 



/i + /i + e and r + — > e + e + p , but no others. In particu- 
lar, /i — > ej is forbidden. Using 

2 It 1, 2\2 

B{t + -> e ) = 4 4 B(r -» /jw), 

the experimental up per limit of 2.3 x 10~ 8 yields the 
bound [Hi mjmj/^+m^ > 22 GeV (174 GeV/v) 
on the masses of ip® 2 = (4>i i 02 ) A/2- 

Since the Higgs singlets x an d ?7 which support the 
neutrino tribimaximal mixing under T7 also transform 
under U(1)b~l, this model can be tested at the LHC by 
discovering the Z' BL (= Z 1 ) gauge boson. The partial 
decay rates of Z 1 to the usual quarks and leptons are 
easily calculated. Let To = g B _ L rnz' /12tt, then T q — 

(6)(3)(i/3) 2 r , r t = (3)(-i) 2 r , r„ = (3)(-i) 2 (i/2)r . 

As for Z 1 — > V'i2 7 / ; 2i; it has the effective partial rate 
~ (2)(-2) 2 sin 4 6»(l/4)r , where sin6» is an effective 
parameter accounting for the mixing of 2 to x an d 
(with the help of a _B — L = singlet Sj ~ 3). Using 
Eq. (18), we find their signature decays to be given by 

i>i,2 t + h~, r-e + , iP h2 -+ r-/i + , r+e", (19) 

resulting in Z' leptonic final states such as r~T~/i + e + for 
example. In addition to being crucial for neutrino tribi- 
maximal mixing to work under T7, the U(\)b-l gauge 
symmetry is seen to provide also the means of verifying 
its predicted interactions. If the singlet neutrinos vf are 
light enough, they can also be produced by Z' decay as 
discussed in Ref. [lij]. The mass eigenstates of v£ are 
given by Eq. (13). Their decays into 4>\^ and leptons, 
and the subsequent decays of ^1,2 to leptons (resulting 
in six leptons in the final state) will then give a complete 
picture of tribimaximal mixing in this model. 

We now study in detail the process qq — >• Z' — >• 
V'lV^ + ip2$i (assuming mj = m^) with the subsequent 
decays ip —> r~e + and ip t~ p + at the LHC with 
E cm — 14 TeV. We consider only the leptonic decay 
modes of the r~, with branching fraction 17.4% to ei- 
ther e~ or \i~ . The collider signature of such events 
is e + [i + l~l~ plus missing energy, where I = e, fi. The 
dominant backgrounds yielding the same signature arc 
the processes (generated by MadEvent/MadGraph [l8l]): 

wwz ■. PP -> w+w-z, W ± -> e^u, Z -> t + r, 

ZZ :pp^r ZZ, Z -> £+£~ , Z -> t+t" , -> £* w, 
ti:pp^ti^ &(-> £")&(-)• £+)TV + VLV PL* -> ^1/, 
Zto : -> Z6(^ £")&(->• £+), Z £+£", (20) 

where £ — e, p. Other SM backgrounds, e.g. ZZZ and 
WWWW , occur at a negligible rate after kinematic cuts, 
and are not shown here. We require no jet tagging and 
consider only events with both e + and p + in the final 
state. The first two processes are the irreducible back- 
ground, while the last two are reducible as they only 
contribute when some tagged particles escape detection, 
carrying away small transverse momentum (pr) or falling 
out of the detector rapidity coverage. 



TABLE I: Signal and background cross sections (fb) before 
and after cuts for four (m Z ',rn^) (GeV) benchmark points: 
(A) (1000,100), (B) (1500,100), (C) (1000,300), and (D) 
(1500,300). The "no cut" rates correspond to all leptonic 
decay modes of r~ after e + /i + identification, the "basic cut" 
and "Ht cut" rates are obtained after imposing Eq. (|21[) and 
Eq. Q22p . respectively, whereas the u x Ti > 0" rates are ob- 
tained after the r~ reconstruction cuts. The bottom row 
shows the cut acceptance (Acut)- 





(A) (B) (C) (D) 


ti WWZ ZZ Zbb 


no cut 
basic cut 
Ht cut 

It > 


5.14 0.98 2.57 0.72 
1.46 0.066 1.05 0.36 
1.41 0.065 1.04 0.36 
0.69 0.032 0.52 0.18 


1.22 0.21 27.11 2.99 
0.16 0.02 0.0052 0.024 
0.08 0.006 0.0 0.0 
0.015 0.002 0.0 0.0 


■A-cut 


13.4% 3.2% 20% 25% 


1.2% 1% 0.0 0.0 



Our benchmark points are chosen as follows: my = 
1000 (1500) GeV, = 100 (300) GeV, Qb—l = 9 = 
e/ sin9w, and sin 2 8 = 0.2. In our analysis all events are 
required to pass the following basic acceptance cuts: 

Pr'f > 50 GeV, p^f > 20 GeV, \r)t\ < 2.5, 
AR U , > 0.4, E T > 30 GeV, (21) 

where (1-4) in the superscript index is the px order of the 
charged leptons. Ai?^ is the separation in the azimuthal 
angle (4>) - pseudorapidity (rj) plane between i and j, 

defined as ARij = y (rji — ijj) 2 + (<fri — 4>j) 2 ■ We also 
model detector resolution effects by smearing the final- 
state energy. To further suppress the SM backgrounds, 
we demand 

Ht = J^Pt, »+ Et > 300 GeV, (22) 

i 

where i denotes the visible particles. Figure Ufa) shows 
the normalized Ht distribution of both signal and back- 
ground before the Ht cut. The signal spectrum exhibits 
an endpoint around the mass of Z', about 1 TeV, but 
with a long tail due to the Z' width and detector smear- 
ing effects. Table Q] displays the signal and background 
cross sections (fb) before and after cuts. The cut ac- 
ceptance (Acut) increases with m^p as heavy scalar decay 
generates hard leptons and large Ht- For a light ip and a 
heavy Z' (e.g. the benchmark B), A cu t decreases as the 
two charged leptons from the light scalar decay are very 
much parallel and fail the AR separation cuts. 

To reconstruct the scalar ip, we adopt the collinear ap- 
proximation that the charged lepton and neutrinos from 
t decays are parallel due to the large boost of the r. Such 
a condition is satisfied to an excellent degree because the 
r leptons originate from a heavy scalar decay in the sig- 
nal event. Denoting by x Ti the fraction of the parent r 
energy which each observable decay particle carries, the 
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FIG. 1: (a) Normalized distribution of Ht\ (b) Distribution 
of the invariant mass of the e + and reconstructed r~pair; (c) 
Distribution of the mass of the reconstructed Z'; (d) The 5cr 
significance contours in the plane of mz' and sin 2 6. 



transverse momentum vectors are related by [19( 

P T = (l/x T1 -l)p 1 + (l/x T2 -l)p 2 . (23) 

When the decay products are not back-to-back, Eq. (|2"3")l 
gives two conditions for x Ti with the r momenta as p\ /x Tl 
and p2/x T2 , respectively. We further require the calcu- 
lated x Ti to be positive to remove the unphysical so- 



lutions. There are two possible combinations of e + £~ 
clusters for reconstructing the scalar ip and gauge boson 
Z'. To choose the correct combination, we require the 
e + £~~ pairing to be such that AR e + e ~ is minimized. The 
mass spectra of the reconstructed ip and Z' are plotted 
in Fig. |Tfb) and (c) , respectively, which clearly display 
sharp peaks around my, and mz>- In Fig. [IJd), we show 
the 5a discovery contours in the plane of mz 1 and sin 2 9 
by requiring 8.5 (5) signal events for an integration lumi- 
nosity of 100 (10) fb _1 respectively. The regions above 
those curves are good for discovery. 

In the quark sector, if we use Qi — (u, d)i ~ 3 and 
u i > di ~ la i — 1)2,3 as we assume for the charged lep- 
tons, we again obtain arbitrary quark masses, but no 
mixing. To have realistic mixing angles, the residual Z% 
symmetry has to be broken. 

Non-Abelian discrete symmetries have been success- 
ful in explaining the tribimaximal mixing of neutrinos, 
but not their masses. However, they are very difficult 
to test experimentally. In this paper, by combining T-j 
and U(1)b-l, we show how a simple renormalizable two- 
parameter neutrino mass model of tribimaximal mixing 
can be constructed, with verifiable experimental predic- 
tions. The key is the possible discovery of Z 1 at the TeV 
scale, which then decays into neutral Higgs scalars, whose 
subsequent exclusive decays into charged leptons have a 
distinct flavor pattern which may be observable at the 
LHC. 
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